Femtosecond filament at ion in air and higher-order nonlinear it ies 
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According to a recent experiment, the instantaneous electronic Kerr effect in air exhibits a strong intensity 
dependence, the nonlinear refractive index switching sign and crossing over from a self-focusing to a de-focusing 
nonlincarity. A subsequent theoretical work has demonstrated that this has paradigm-changing consequences 
for the understanding of filamcntation in air, so it is important to subject the idea of higher-order nonlinear- 
ities to stringent tests. Here we use numerical modeling to propose an experiment capable of discriminating 
between the standard and the new intensity-dependent Kerr-effect models. © 2010 Optical Society of America 
OCIS codes: 320.2250, 320.7110. 



Over the last sesquidecadc research into femtosecond 
filamcntation in air has produced a widely accepted pic- 
ture in which a dilute plasma, generated in the high- 
intensity optical field, is a major determinant of the light 
pulse dynamics. For gaseous media, this scenario is based 
on the notion of a dynamical balance between the Kcrr- 
induced self-focusing and the counteracting de-focusing 
caused by the free electrons [T]. Higher-order nonlinear- 
ities have previously been considered before as an addi- 
tional mechanism regularizing the self-focusing collapse. 
In particular, phenomenological Kerr-effect saturation at 
high intensities (see e.g. [2H1] and Sec 3.4. of [5]) was ex- 
plored in filament modeling. In several papers quantum- 
mechanical simulations of the electronic polarization re- 
sponse (e.g. [5] and references therein) were used to study 
the interplay of the saturated Kerr and plasma contribu- 
tions. Nevertheless, in the absence of experimental evi- 
dence for these effects, the consensus in the filamcntation 
community was that the standard self-focusing Kerr ef- 
fect was sufficient to describe the electronic nonlinearity. 
Moreover, many experiments clearly showed the pres- 
ence of plasma in filaments, and this was accepted as 
the mechanism clamping the intensity [7J. 

Recent experimental results by Loriot et al. [HUHj have 
the potential to change this long-standing paradigm. 
They found that the instantaneous Kerr nonlinearity ex- 
hibits a strong intensity dependence, and that it not 
only saturates, but even changes from self-focusing to 
de-focusing at large enough intensity. As a consequence, 
the plasma is not necessary for the arrest of the self- 
focusing collapse, because the Kerr-induced de-focusing 
at high intensity can alone be sufficient. Inspired by these 
findings, the recent modeling work of Ref. [10] proposed 
that contrary to our current belief, ionization-free fila- 
mentation is possible in gases. 

This has been a surprising development to say the 
least, and one which poses many questions. The first 
one is of course, why the earlier experiments did not 
detect such a strong intensity dependence of the Kerr 
effect? Further, is it just an accident that the cross-over 
from self-focusing to de-focusing occurs in the vicinity 
of the intensity range where plasma starts to form, or is 



there a connection between the two? In particular, cen- 
tral to the interpretation of the Loriot's experiment is 
the assumption that the optical response of the plasma 
is isotropic, which is certainly true for the simple Drude 
model. However, ultrafast birefringence is know to occur 
for free electrons in semiconductors due to anisotropic 
state filling [11] , and if a similar effect arises for the free- 
carrier plasma in air it could change the interpretation 
of the experiment. And perhaps the most fundamental 
question is what is the microscopic origin of the sign- 
change in the Kerr effect if not due to plasma? 

Since the observations of Loriot et al. have such large 
implications for the physical picture of filamentation in 
gases it is important that they be corroborated in simi- 
lar experiments, and that the occurrence of higher-order 
nonlincarities be subjected to independent experimental 
testing. This may not be easy, because accurate quantita- 
tive measurements in filaments are notoriously difficult. 
It is also conceivable that the Kerr intensity-dependence 
could be partially masked by the plasma induced de- 
focusing. In fact, this may very well be the reason why 
this effect has not been observed earlier. It is therefore 
the view of the present authors that this problem needs 
to be addressed with an utmost caution. Before making 
conclusions with regards to the filamcntation dynamics, 
ways should be found to independently test the interpre- 
tation of the Loriot's experiment. As a first step in this 
direction, based on numerical simulations we propose 
the design of an experiment which should be capable 
to discriminate between the two model candidates: i.e. 
between the standard Kerr, and the higher-order Kerr 
effects. 

The proposed experiment is simple in principle. One 
should use a modest pulse energy to generate third- and 
fifth-harmonic radiation, and measure their respective 
conversion efficiencies. To make this feasible, the funda- 
mental frequency should be "shifted" to a longer wave- 
length, perhaps by using an OPA, such that the third and 
fifth harmonic energies can be simultaneously measured. 
In order to minimise the propagation effects, the focusing 
geometry should be tight, resulting in a short interaction 
volume. The aim is to compare the fifth-harmonic gen- 
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eration to that of the third harmonic, and decide if the 
result is in keeping with the predictions of either of the 
competing models. The rationale behind this is the sim- 
ple fact that the higher-order Kerr effect should result 
in a much stronger harmonic generation. Indeed, in the 
standard model only the third harmonic is produced di- 
rectly from the fundamental, while the fifth (and higher) 
harmonics arise via a cascade from the third. On the 
other hand, the fifth harmonic can be produced directly 
from the fundamental in the higher-order Kerr model. 

In our numerical experiments we use the Unidirec- 
tional Pulse Propagation Equation simulator [12] which 
works directly with the real electric field (as opposed to 
amplitude envelopes) , and therefore naturally and with- 
out any modifications captures all frequency-mixing in- 
teractions. We assume the 50 fs input pulse with a central 
wavelength of 1330nm and specify its nominal intensity 
and beam waist which would be achieved in a linear fo- 
cus. The results presented in this paper were obtained for 
the nominal focused beam waist of 50 microns, and the 
nominal in-focus intensity between 10 17 and 10 18 W/m 2 . 
This tight focusing results in the interaction region only 
several centimeters long and is meant to suppress the 
possible propagation- and plasma-related effects. For the 
plasma generation, we use the so-called PPT ionization 
rates [T3] parametrized as in Ref. [T3J . Besides the instan- 
taneous Kerr effect, we include the usual delayed Raman 
component [12) with the instantaneous-to-delayed ratio 
of one half. 

It should be emphasized that the only difference be- 
tween the comparative simulation runs is the different 
inputs specifying the nonlinear index as a function of the 
instantaneous value of the electric field. For the higher- 
order Kerr effect model, we have (with suitable electric 
field units) An = n 2 E 2 + n 4 E 4 + n 6 E 6 + n 8 E s with the 
7&2fc values for air fixed according to [9]. The standard 
model shares the same n 2 and has n 4 = n§ = ng = 0. 

Let us first look at the typical spectra generated by the 
two models. Figure [T] shows the results for the nominal 
focal intensity of 10 18 W/m 2 (the actual intensity reached 
in the interaction volume is somewhat smaller). As ex- 
pected, we can see a drastic difference. First, in the new 
model the third-harmonic generation is almost two or- 
ders of magnitude stronger. Second, the fifth harmonic 
generation is four to six orders of magnitude stronger. 
Thus, if the new model is correct the fifth harmonic 
should be easily detected. 

To quantify the above observation further, in Fig. [5] 
we show how the third- and fifth-harmonic conversion 
efficiencies evolve with the increasing pulse energy (or, 
cquivalcntly, with the nominal focal intensity) . The stan- 
dard model produces the 3H efficiency approaching frac- 
tion of a percent which is in agreement with both exper- 
iments (see e.g [15l[T6]) and previous simulations [17] . 
On the other hand, the 3H efficiency in the higher-order 
model can be as high as several percent! To our best 
knowledge this has not been observed yet. 

However, one could argue that a mere order of mag- 
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Fig. 1. Angle-integrated spectrum showing third- and 
fifth-harmonic generation. The standard and higher- 
order Kerr effect models are compared for a tightly 
focused femtosecond pulse. The new model predicts 
stronger third-harmonic generation and several orders 
of magnitude more powerful fifth-harmonic radiation. 
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Fig. 2. The harmonic conversion efficiency as a function 
of nominal focal intensity /energy. Full and open sym- 
bols represent the third and fifth harmonics, and squares 
and diamonds represent the higher-order and standard 
model, respectively. 

nitudc difference in the third-harmonic conversion effi- 
ciency might be difficult to detect experimentally, be- 
cause any real experiment will only probe "a single model 
version." That is why we now turn to the fifth harmonic 
conversion efficiency Figure [5] shows that here the dif- 
ference reaches five orders of magnitude, and this should 
indeed be possible to identify in the experiment. In fact, 
if the standard model is correct, the fifth harmonic might 
be hardly detectable as it should be seven orders of mag- 
nitude weaker than the fundamental. If on the other 
hand the new model is correct, the fifth harmonic should 
be easily detected on the level of up to a few percent of 
the fundamental energy. 

Another potentially useful telltale sign which an ex- 
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Fig. 3. Ratio of the fifth- to the third-harmonic energy 
generated in the tightly focused femtosecond pulse. The 
standard model ratio roughly corresponds to an expo- 
nential decrease between the harmonic orders. In con- 
trast, the new model predicts that the energy carried by 
the fifth harmonic can be a mere order of magnitude less 
or even comparable to that in the third harmonic. 



periment may look at is the ratio between the fifth and 
third harmonic production. This is shown in Fig[3l While 
the standard model exhibits a small ratio reflecting the 
fact that the process is a cascade and the fifth harmonic 
production feeds on the third harmonic energy, the two 
harmonics become almost comparable in strength in the 
new model. This occurs as soon as the focal intensities 
reach values beyond the turn-over point of the higher- 
order Kerr effect. In fact, it is easy to envision that once 
the electric field amplitude increases beyond the turn- 
over point, the nonlinear index change will exhibit two 
peaks within each half-cycle of the driving field. This in 
turn means that the nonlinear index spectrum contains 
roughly the same amounts of second and fourth harmon- 
ics giving rise to comparable strengths of the third and 
fifth harmonic components in the nonlinear polarization. 
In a nutshell, this is the reason why the new higher- 
order Kerr model should produce extremely strong fifth- 
harmonic radiation. 

Let us summarize our findings. Already the present 
calculations suggest that the standard model is in a bet- 
ter agreement with previous experiments on the third- 
harmonic generation. It exhibits the conversion efficien- 
cies comparable with the measured values. In contrast, 
the new model predicts significantly stronger third har- 
monic, with its energy reaching up to several percent of 
the fundamental. To our knowledge such a strong third- 
harmonic generation has not been observed. However, we 
should keep in mind that due to the model parameter un- 
certainties the possibility can not be excluded that the 
simulated efficiency value may suffer from a systematic 
deviation. 

That is why we propose an experiment which should 
compare the fifth and third harmonic generation. We 



have shown that if the higher-order Kerr effect model 
is correct, the experiment will detect significant fifth- 
harmonic generation. Moreover, at higher pulse energies 
the relative strength of the fifth to the third harmonic 
should reach values of the order of 10 _1 . On the other 
hand, if the standard model is correct, this ratio should 
be about four to five orders smaller. This is a signifi- 
cant qualitative difference, and it should be obvious from 
the experimental results which of the competing models 
should be preferred. 

In our view, the importance of the independent cor- 
roboration of the experiment by Loriot et al. can not be 
overemphasized. Their results imply literally a change of 
paradigm in our understanding of femtosecond filamen- 
tation, put in doubt concepts such as plasma-mediated 
intensity clamping, and essentially invalidate the very 
model of femtosecond filamentation that has been used 
for over a decade. It is thus our aim and hope that this 
work will motivate and inspire experimental efforts to 
solve this exciting and important problem. 

This work was supported by Air Force Office of Scien- 
tific Research under contract FA9550-07-1-0010. 
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